Master equation for a particle coupled to a two-level reservoir 
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Abstract 



We study the quantum dissipative dynamics of a particle coupled linearly 
to a set of two-level systems (the heat bath) via the master equation method 
which we extract from the path integral formalism independently from the 
form of the bath spectral density. We compare our results with the standard 
models based on bosonic heat baths showing their main differences and sim- 
ilarities. In particular, we study special forms for the spectral density of the 
bath which give results quite different from the standard models. 
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The study of quantum open systems is one of the most important issues of statistical 
quantum mechanics because of its application to many different subjects such as thermaliza- 
tion (tendency to thermal equilibrium), lost of quantum coherence and quantum brownian 
motion. This kind of theoretical approach has been used with success in order to study 
macroscopic quantum tunneling [TJ,|2|], evolution of the early universe |3| and measurement 
theory ||, among others. 

The field of quantum brownian motion had a major revival after the work of Caldeira 
and Leggett [0 on the study of a particle coupled to a set of harmonic oscillators via the 
method devised by Feynman and Vernon || using path integrals. The path integral method 
allows the trace over the oscillators variables (the heat bath) and the study of the coordinate 
of interest. As it is already well known the dissipation and diffusion of a particle coupled 
linearly to a heat bath depend essentially on the spectral density, J(u), of the heat bath 
(which is related to the response function of the reservoir at some frequency uS). In the 
Caldeira-Leggett model [|1| the dissipation is supposed to be ohmic, that is, the spectral 
density of the bath has the form, 

J (u) = i] u 9(uj c -u), (1) 

where u c is a suitable cut-off frequency. With this special choice for the spectral density 
Caldeira and Leggett showed that the semiclassical dynamics of the system is described by 
the Langevin equation with a dissipation coefficient r/o- 

However, the study of open quantum systems has a long history with many branches. 
One of the main streams of research is based on the study of master equations (such as the 
Fokker-Planck equation) for the dynamical evolution of quantum systems which is used in 
quantum optics |||7[]. In the seventies, Lindblad gave a general mathematical criteria for 
the classification and existence of such master equations ||. More recently, the methods of 
master equations has been used extensively in the literature P JT0|JT^] . Actually, Caldeira 
and Leggett obtained the equation for the evolution of the reduced density matrix, p s , (the 
density matrix after the trace over the oscillators bath) with the ohmic spectral density 
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(1) and in the high temperature limit (ksT ^> hu c ) which, in the position representation 
(p s (x,y,t) = (x\p a (t)\y)), reads II, 



dp s {x,y,t) 

in 

dt 



— (—-—) -ihYix-y) 

2m \ dx 2 dy 2 1 \dx dy 



+ V R (x) - V R (y) - % D pp (x - yf 



Ps(x,y,t), 



(2) 



where m is the mass of the particle, Vr is the renormalized external potential and 

2rj k B T 



D 



pp 



r _ Vo 

m 



are the diffusion coefficient in momentum space and the effective dissipation coefficient, 
respectively. 

In this paper we address the problem of the master equation for a different kind of bath 
which was proposed recently in the literature by Caldeira et al. [I3| and has distinct features 
from the models described above. In the model proposed in Ref. |l3j the particle interacts 
with a heat bath composed of two-level systems. The hamiltonian of interest can be written 
as a sum of three terms, 



H = H s + Hi + H b 



where H s is the Hamiltonian for a particle in an external potential V(x) 

-,2 



H* 



Hi gives the interaction between the particle and the bath, 



Hi 



X O r 



and the bath hamiltonian is, 



H b 



E?lUJ n 



where a Xn and a Zn are the usual Pauli matrices and J n is the coupling constant for the 
particle with the n th two-level system. 



As it was shown in ]T3|], the dynamics of the particle can be evaluated in the path integral 
formalism via the superpropagator, J, |J. Assuming that initially at some time t = the 
particle and the bath are decoupled the time evolution of the system can be written as, 



Pa{x,y,t) 



where the superpropagator is given by, 



dx / dy J(x,y,t | x ,y ,0) p s (x ,y ,0) 



(3) 



J(t,0) 



Vx(s) I Vy(s) exp — 

xq Jyo h 



S[x(s))-S[y(s)) 



F[x(s),y(s)} 



(4) 



In the last equation, F [x(s),y(s]) is the influence functional ||, and S [x(s)] is the free 
action for the particle, 

mx 2 (s) 



S[x(s)} 



ds 



~ V[x{s)\ 



In this paper we will be interested in the case where the external potential is harmonic, that 
is, V [x(s)] 



^mfij x 2 (s). 



In the weak coupling limit the influence functional reads JTB 



F[q(s),Q(s)} = exp -- 



x exp 



ft rs\ 

/ dsi / ds 2 2 q(si) Q(s 2 ) v( s i ~ s 2) 
Jo Jo 

ft fSl 

/ dsi / ds 2 q(si) q(s 2 ) v(s\ - s 2 
Jo Jo 



(5) 



where, 



q{s) 



x ( s > - y( s ) 



Q(s) = \ W) + y(s)) 



(6) 



and the kernels in (5) are given in terms of the spectral density as, 



VKS) 



du J{oj) cos us 



(7) 



and 



POO 

/ du JUS) tanh 
Jo 



2k B T 

For this particular model the spectral density is written as, 



sin us. 



(8) 



n ^ 

Since the double path integral in (4) is quadratic it can be evaluated completely in terms 
of the initial and final coordinates given by the transformation (6), the result is : 



J(t, 0) = Z(t) exp - 



(hi Qq + b 2 Q) q - (6 3 Qo + h Q) Qo 



1 

x exp — — 
n 



an q 2 + a\2 Qo <? + O22 Qo 



(9) 



where the coefficients b n (n — 1, 4) and a nrn (n,m = 1,2) are given by the solution of the 
equations : 

Cp K ( s) f s 

with the boundary conditions, 

Ki(s = 0) = K 2 (s = t) =1 Ki(s = t) = K 2 (s = 0) =0. 



Then the unknown functions in (9) are defined as, 



Ht) = m ds 

b 3 {t) — M — — 
as 



s=t 



s=0 



b 2 (t) = M 
b 4 (t) = M 



dn 2 

ds 
dn 2 
ds 



s=t 



and, 



1 ,-t ,-t 

a nm (t) = —— £ — / ds 1 / ds 2 K n {si) v{s x - s 2 ) K m (s 2 ). 
1 + d nm Jo Jo 

Also the normalization factor in (9) is given by, 



Z(t) = 



2nh 



Using the above equations we can derive a set of identities, 



h a 12 



1 

M 



(10) 



bi b 2 b A 

(the dot over the functions denotes a derivative with respect to t), which we shall use latter 
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The master equation is obtained using the simplified method proposed by Paz [IT . 
Taking the derivative of (3) with respect to t it is straightforward to show that, 



dp a (q,Q,t) 
dt 



Z * i. r\ 1 • 

77 + ^2?Q - rang 



i . r+00 
+ - bi q / dg 

/I J— CO 7 —CO 

* L ^ 1 • 

- 64 Q + - a i2 <? 
T h / d<?o 

/2 J —00 



h J -00 



Ps(q,Q,t) 
dQo Qo J(t,0) p s (qo,Qo,0) 

+00 /"+oo 

d<? / dQo qo J(t,0) p s (qo,Qo,0) (11) 

dQo go Qo J(t,o) p s (qo,Qo,o) 

dQo ql J(t,0) p s (qo,Qo,0) 



Although (11) seems quite complicated we can simplify it by taking derivatives of (9) 
with respect to the end points, q and Q, which lead to the following relations, 



Qo J(t, 0) 



b 2 1 d 

t q + lh r ^tt 

o 4 o 4 oQ 



J(t,0) 



(12) 
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J(t,0). 



By substituting (12) in (11) and using (10) and (3) we find after some straightforward 
algebra that, 

h 2 ( d 2 d 2 



dp B (x,y,t) 

in : 

dt 



2M \dx 2 dy 2 



+ 



MQ 2 R 2 



x-y 



Pa(x,y,t) 



( d d\ 

- in T(t) (x - y) I — - — I p s (x, y, t) 

- i Dpp(t) (x - y) 2 p s (x,y,t) 

+ h [D xp (t) + D px (t)] (x - y) l-^- + p s (x,y,t) 



(13) 



is the master equation for the reduced density matrix of a particle coupled to a two-level 
reservoir. Observe that the term which represents the diffusion in real space, D xx , is absent 
in this model since the kernels in (7) and (8) depend only on the time interval Si — s 2 (see 
Eq. (5)). This is a general feature of stationary kernels flTOf . 
The coefficients that appear in (13) are given by, 
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is the renormalized frequency for the harmonic potential, 

bi 



r(t) 



2 61 



(14) 



is the dissipation coefficient, 



D pp (t) 



4 an b 2 + a 12 



b_i 



2 an + ai2 



(15) 



is the diffusion coefficient in momentum space (it gives the fluctuation in x 2 ) and 



D xp (t) = D px (t) 



Ql2 

64 



h Qi2 
61 64 



4 an 



(16) 



is the diffusion coefficient which mix the real space and the momentum space (it gives the 
fluctuation in xp + px) . 

The first observation is that the kernels in (7) and (8) do not have the same form as for 
the case of the oscillator bath |l| unless we choose a temperature dependent spectral density 
of the form, 



J\(ui) = rji to coth 



2k B T 



(17) 



as it was first realized in Ref. jnj (it reduces to the spectral density (1) when T — * 0). With 
this choice the bath of oscillators is indistinguishable from the bath of two-level systems. 
For example, in the high temperature limit (ksT ^> hu c ), using (17) we find, 



2r] 1 k B T 



5(s) 



7]{s) 



dS(s) 



h ds 
In this case the diffusion and the dissipation are markovian while at low temperatures the 
diffusion is clearly non-mar kovian. Moreover, the coefficients that appear in (14) — (16) are 
given by, 



T(t) 
as expected. 



ni 

m 



D xp {t) = D px {t) =0 



D pp {t) 



2r] 1 k B T 

h ' 



(18) 



However, if we insist in using the ohmic spectral density, Eq.(l), the kernels in (12) and 
(13) assume, at high temperatures {k B T ^> huj c ), a non-local form, namely, 

/ \ 1 7] h l 

and the dynamics is non-markovian. It is interesting to notice that in the low temperature 
limit {k B T <C hu c ) the use of this ohmic spectral density gives the following result for the 
kernels : 

1 d5(s) 



s 2 1 ds 

and the only time independent coefficient in the master equation is the dissipation one : 
F(t) = This shows that the two- level system and the oscillator bath cannot be naively 
mapped into each other with this spectral density. 

Another interesting possibility is a flat spectral density of the form, 

J 2 (u) = r] 2 6(Lu c -Lu), (19) 

which is very common in applications in quantum optics ||[14| (it should be noticed that 
rj2 is dimensionally different from r/ , 771). In this case the diffusion is purely markovian 
(without memory) at any temperature while the dissipation coefficient presents a transition 
from markovian to non-markovian depending on the temperature scale, that is, substituting 
(19) in (7) and (8) and taking the high temperature limit, we find, 

r] 2 h dS(s) 

and the coefficients in the master equation are given by, 

r(t) = D xp (t) = D px (t) =0 D pp (t) = V2 . 

Observe that in this case the dissipation coefficient has temperature dependence (it decreases 
with the temperature) and the diffusion coefficient is temperature independent. This result 
must be compared with (18) where the dissipation does not depend on the temperature but 
the diffusion increases with the temperature. 
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In conclusion, in this paper we obtained the most general master equation for a har- 
monic oscillator coupled to a reservoir consisting on two-level systems. We show that the 
kernels which generate the diffusion and dissipation coefficients in the master equation have 
a different form from the standard cases of heat baths composed of harmonic oscillators. We 
confirm the results of Caldeira et al. [I3"| that the results obtained with an oscillator bath 
can be recovered from a two-level system bath by a temperature dependent spectral density 
(which reduces to the ohmic one when T — * 0). We showed that the usual ohmic spectral 
density produces a different kind of behavior from the usual oscillator model [[|]. We also 
studied the case of a flat spectral density (usually used in quantum optics) which produces 
a different kind of diffusion and dissipation for the quantum dynamics of the particle. 
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